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We consider the problem of distinguishing with minimum probability of error two optical beam 
splitter channels with unequal complex-valued reflectivities using general quantum probe states 
entangled over M signal & M' idler mode pairs of which the signal modes are bounced off the 
beam splitter while the idler modes are retained losslessly. We obtain a lower bound on the output 
state fidelity valid for any pure input state. We define Number-Diagonal-Signal (NDS) States to 
be input states whose density operator in the signal modes is diagonal in the multimode number 
basis. For such input states, we derive series formulae for the optimal error probability, the output 
state fidelity, and the Chernoff-type upper bounds on the error probability. For the special cases 
of quantum reading of a classical digital memory and target detection (for which the reflectivities 
are real- valued) , we show that for a given input signal photon probability distribution, the fidelity 
is minimized by the NDS states with that distribution and that for a given average total signal 
energy Ns, the fidelity is minimized by any multimode Fock state with Ns total signal photons. 
For reading of an ideal memory, it is shown that Fock state inputs minimize the Chernoff bound. 
For target detection under high loss conditions, a no-go result showing the lack of appreciable 
quantum advantage over coherent state transmitters is derived. A comparison of the error probability 
performance for quantum reading of number state and two-mode squeezed vacuum state (or EPR 
state) transmitters relative to coherent state transmitters is presented for various values of the 
refiectances. While the nonclassical states in general perform better than the coherent state, the 
quantitative performance gains differ depending on the values of the reflectances. The experimental 
outlook for realizing nonclassical gains from number state transmitters with current technology at 
moderate to high values of the refiectances is argued to be good. 

PACS numbers: 42.50.Ex, 03.67.Hk 



I. INTRODUCTION 

Consider, for b € {0, 1}, an optical beam splitter channel 
taking an input signal mode annihilation operator a\„ to 
the output mode annihilation operator Cout in the manner 
of Fig. 1 under the mode transformation;- 




Here, 7'b = ^/Rh and it = \/7b are real and non-negative 
field reflectivities with i?b and the corresponding re- 
fiectances and transmittances with i?b + 71, = 1. We 
will assume Rq < Ri throughout the paper. We have 
included, for later purposes, the input and output anni- 
hilation operators e\n and Cout of the other (environment) 
mode incident at the beam splitter, whose input is as- 
sumed to be in the vacuum state. Denoting the quantum 
channels induced on the signal mode by the above beam 
splitter transformation by fbi b = 0/1, and assuming that 
each of these channels has equal a priori probability, we 
consider the following general strategy to discriminate 
these two channels with minimum error probability illus- 
trated in Fig. 1:- A quantum state source S prepares a 
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pure state on the system consisting of M 'signal' opti- 
cal modes {ar^}^=i and M' 'idler' modes {C'j^Ci, al- 
lowing for any entanglement across the signal modes and 
between the signal and idler modes. The M signal modes 
are passed through £h with the idler modes unchanged, 
giving rise, in general, to mixed density operators 

p,^£f' ^\d^'''{\tP){^\) (2) 

on the signal-idler Hilbert space depending on the value 
of b - here, id denotes the identity channel on the idler 
modes. This joint state is now measured at the detec- 
tor D using the Helstrom quantum measurement that 
yields minimum error probability in distinguishing the 
two states [T] . In view of the fact that loss is ubiquitous 
in quantum state transmission, processing and detection, 
and because of the well-known sensitivity of nonclassical 
states to loss, we will in general assume that r^, < 1. 
On the other hand, at room temperature and at opti- 
cal wavelengths ^ 1 /xm, the average number of ther- 
mal noise photons per space-time mode is ideally around 
10~^^, allowing us realistically to neglect it, although it is 
known that interesting nonclassical gains are obtained in 
high-loss high-noise conditions such as those in 'quantum 
illumination' \2^. 

Many interesting problems are encompassed by the 
above model. The case of rg = (and possibly n <C 1) 
corresponds to a target detection scenario (identical to 
quantum illumination [2 but for the absence of ther- 
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FIG. 1: Schematic of setup for determining which of two beam 
sphtter channels (indexed by b) of the form of ([T| is present 
within a black box B. A quantum state source S produces a 
pure state of M signal modes {a™}m=i and M' idler modes 
{pin }m'=i of which One signal-idler pair is shown. Each signal 
mode reflects off the beam splitter while the idler mode is 
unaffected. The environment modes {e™}^^i coupling to 
the beam splitter are all in the vacuum state. The detection 
module D performs the minimum error probability (Helstrom) 
measurement for distinguishing the two possible (M + M')- 
mode quantum states corresponding to b = 0/1. 

mal noise) in which b = 0, 1 represents the absence or 
presence of a reflecting target of effective field reflec- 
tivity ri. The more general case of non-zero rg but 
A = Ox — 9q — corresponds to the recently proposed 
quantum reading of a classical digital memory [3]. The 
case of rg — ri < 1,A 7^ models the lossy discrimi- 
nation of channels differing only in phase shift (or path 
length) , and is related to the problem of estimation of a 
continuous optical phase parameter. 

Given this background, the following question may now 
be posed: Among all \ip)'s with average total energy Ng 
in the signal modes, which state minimizes the average 
error probability in distinguishing the channels, and how 
much improvement over a classical state of the same en- 
ergy is attainable? A multimode classical state is a state 
with a non-negative P-representation [HH]. Thus p given 
by 

p = y P{a)\a){a\(fa (3) 

is classical if P{a) > for jo;) an M-mode coherent state, 
so that p is a probabilistic mixture of product coherent 
states. Comparison of the performance of a proposed 
input state to this class of states is important since, as 
their definition implies, the latter are readily prepared by 
modulating laser fields with classical random numbers. 
In particular, we will compare performance with the pure 
coherent state lO)^"''' 0. 

A related version of the general problem of discrim- 
inating two beam splitter channels has been addressed 
recently in |7|. The scenario of [7] differs from ours in 



that it is assumed, referring to Fig. 1, that the com- 
bined state of Cin and din may be chosen subject to a 
total energy constraint and that both the beam splitter 
output modes dout and eout are available for a Helstrom 
detector to measure, in addition to the idler mode b\n- 
As such, the problem of [7] is one of discrimating two 
unitary transformations as opposed to the non-unitary 
channel discrimination problem considered here. Never- 
theless, it is clear that, for M = 1, the minimum error 
probability attainable at given Ng in the setting of [7 is 
not greater than that for our more restrictive setting. On 
the other hand, the case of M > 1 was not considered 
in [7]. We note that the problem in the form addressed 
here is of more relevance to long-range scenarios where 
the black box B of Fig. 1, and therefore, the input and 
output environment modes, are not directly accessible to 
the user who controls the state source S and the detector 
D. It also models quantum reading scenarios where the 
memory itself (e.g., a CD or DVD) need not be modified 
in a major way. When such modifications to the memory 
are contemplated, the scenario of [7| becomes interesting. 

This paper provides several partial answers to the 
broad question posed above. In Section II, we derive a 
lower bound on the Jozsa-Uhlmann fidelity [SI IS] between 
Pq and pi 

J" = Tr^ ^J^/Jkipu/Po (4) 

that is valid for any input state [10] . This further 
yields a lower bound on the minimum error probability 
of distinguishing the output states. In Section III, we spe- 
cialize to input states whose reduced density operator on 
the signal modes is a mixture of multimode number states 
- such states will be referred to as Number-Diagonal- 
Signal (NDS) States in this paper. It was shown recently 
in [TT] that NDS states are optimum input states ac- 
cording to many possible performance criteria for a large 
class of image sensing problems of which the minimum 
error probability discrimination between two beam split- 
ter channels that is the subject of the present paper is a 
special case. We mention that, in the contexts of commu- 
nication and key distribution between two users, a class 
of states related to the NDS states was introduced in 
[Hj and referred to as Photon-Number Entangled States 
(PNES). A PNES is a pure state of a single signal-idler 
mode pair that is diagonal in the number state basis of 
both the signal and idler modes. Since the states we con- 
sider do not have to be number-state-diagonal in the idler 
modes, we believe the term "Number-Diagonal-Signal" 
state is more appropriate. However, it is the case that a 
given single-mode NDS state can be related to a corre- 
sponding PNES by a unitary transformation on the idler 
modes without changing the error performance. 

The NDS states include, but are not limited to, number 
states, the two-mode squeezed vacuum (or EPR) states 
used in [IE], the NOON states [T3j, the related |m :: m') 
states of [14], the pair coherent states (PCS) [H], and 
several other states studied in the context of phase esti- 
mation TB^. For this class of states, we show that po and 
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pi of ([2]) are easily diagonalized. The minimum possible 
error probability of discriminating any two states (with 
equal a priori probabilities) is given by the Helstrom for- 
mula [1^ 
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Po - Pi 



(5) 



where || A ||i is the trace-norm which equals, for self- 
adjoint A, the sum of the absolute values of the eigen- 
values of A. For non-commuting mixed states, Pg is no- 
toriously hard to calculate. For input NDS states, how- 
ever, we show in Section III that the eigenvectors of po 
and pi have a mutual inner product structure that per- 
mits the calculation of as a, in general infinite, series. 
The fidelity J4]) also yields the following upper and lower 
bounds on Pe [12] (the first inequality follows from a 
binomial expansion of its RHS):- 

\:F<\{i-Vi^)<Pe<\^- (6) 

Also of interest are the "Chernoff-type" upper bounds:- 



Pe < Iq{s), [0,1], 



(7) 



where Q{s) = Tt[pqp\ and the Chernoff bound Q [18] 
is the best such bound:- 



Pe < 



(8) 



Chernoff bound among all pure input states of given en- 
ergy. In Section IV. A. 3, we state a fidelity-based error 
probability lower bound for target detection with a gen- 
eral input state. In the limit i?i ^ 1 of large loss, we 
derive a no-go result that rules out appreciable quantum 
advantage over coherent states. 

In Sections IV.B-D, we characterize in detail the er- 
ror probability performance of coherent states, number 
states and the EPR states. We show that the perfor- 
mance of all the multimode Fock states of total photon 
number N s are identical, a fact that helps practical imple- 
mentation. For EPR states, we show that the application 
of the results of Section III yields the same analytical re- 
sults as the methods used in [5^. In Section IV.E.l, we 
compare the quantitative performance of coherent, num- 
ber, and EPR states for some typical reflectance values 
and demonstrate nonclassical gains. In Section IV. E. 2, 
we consider the technological feasibility of achieving the 
nonclassical gains in quantum reading. 



II. LOWER BOUND ON OUTPUT STATE 
FIDELITY 



Consider an arbitrary (M' 4- Af)-mode state 
photon number representation:- 



where 



un|m)7|n)s, 



in the 
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min Q(s). 
se[o,il 



(9) 



The Bhattacharyya bound is Q with s = 1/2. This ter- 
minology was introduced in [19] and the bound was ap- 
plied for the first time in [2]. We show here that, for 
NDS input states, and Q{s) are also calculable as in- 
finite series in general, and have closed-form expressions 
in some cases. 

In Section IV, we apply these results to quantum read- 
ing and target detection, i.e., to the cases where A = 0. 
The input state consisting of M signal and idler mode 
pairs from the output of a parametric downconversion 
process was shown in [3] to yield, in some regions of M 
and Ns , surprisingly better error probability in quantum 
reading than that obtainable from any classical state of 
the same energy [20]. In this paper, we first show in 
Section IV. A that, for any given signal photon number 
distribution, the general fidelity lower bound derived in 
Section I is attained by NDS states. Further, the NDS 
states minimizing the fidelity for a given and M are 
the multimode Fock states of total photon number Ng. 
In Section IV.A.l, we lower bound the error probability 
for quantum reading and target detection via the fidelity 
for general input states. In Section IV. A. 2, for the case 
of i?i = 1, we show that the Fock states give the lowest 



where |m)/ = \mi)i^ \mM')i^,, and \n)s = 

\ni) ® ■ ■ ■ ® \nM)sM are respectively the photon num- 
ber states of the M/M'-mode signal/idler. After aug- 
menting the above state with the vacuum e of the M 
environment modes {einlm^i of ([!]), we may write down 
the Schrodinger-picture evolution corresponding to the 
Heisenberg-picture evolution ([T]) of the joint system to 
one of two pure states as:- 



(b)\ 



(11) 



k;<r 



where the amplitude A^^\^ is given by 



A. 



(b) 



M 

n 




(12) 



In (11), k < n is to be understood as component- wise 
inequality. We then have pb = Tr£;[|V'(''^)(V'^''^|] so that 
the IV''''^-') are purifications [ST] of the respective pb- From 
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([TT]), it follows that the overlap O = |(V'(°)|V''^^)P equals 
O = (13) 



M 



k<n m=l 



™ 1 ^ri^^C^ 



(rori)^ 



(14) 



A// 



E^nEn 

n k<nm— 1 



e™-^ (rori)"--'="(ioti)'=' 



M 



E^n n [e""^(^ori+iot 



X]p„e™^(rori+ioti)" 



(15) 
(16) 

(17) 



where l^in)/ is any orthonormal set of states on the idler 
modes - the above equation is essentially a Schmidt de- 
composition [21] of IV')- This is a wide class of states, 
anany of which have been intensively studied in quan- 
tum optics and have interesting applications in quantum 
information and metrology [2 IT^HT^ . 

We now show that, for NDS inputs, po and p\ have a 
form helpful for calculations. After propagation through 
the beam splitter channel we obtain, as in (11), the pu- 
rifications 



k<n 



(22) 



E( E ^n<ll'^n)/|n-k)s)|k)^ (23) 

k n:n>k 

:El^k'>|k)s, (24) 



where p„ = 



is the multimode photon prob- 



ability distribution in the signal modes. In going from 
(161 to ( 17 ) , we have re-ordered the sum of (161 over the 



single index n = X]m=i '^m with p„ being the probability 
distribution of the total signal photon number:- 



Vn 



E P"' 



(18) 



Uhlmann's theorem [SKSIHT] states that the fidelity T is 
the maximum overlap over all purifications of po and pi, 
so that we have the lower bound 



0<T 



(19) 



for the O of (17). For any proposed input state ji/'), the 



overlap 0(|?/'))can be calculated via (17) in terms of the 
signal photon probability distribution oi\i\)). Application 
of the inequalities ( 19 ) and ([6| yields the following lower 
bound on P^:- 



(20) 



The cases of quantum reading and target detection for 
which A = are further developed in Section IV. A. 



III. NUMBER-DIAGONAL-SIGNAL (NDS) 
STATES: MINIMUM ERROR PROBABILITY, 
OUTPUT STATE FIDELITY, AND 
CHERNOFF-TYPE BOUNDS 

We define Number-Diagonal-Signal (NDS) States to be 
states whose reduced density operator in the M signal 
modes is diagonal in the product number (Fock) state 
basis. A pure NDS state then has the representation 



(21) 



where 



(b)N 
k / 



E ^-^^ 



(b) 



On)/ n - 



(25) 



i>k 



are un-normalized joint signal-idler states. Since the 
{|0n)/} arc an orthonormal set by the NDS state as- 
sumption, we see that the orthogonality relations 



(26) 



hold because the and |V'k"0 of ( |25[ ) are termwise 

orthogonal when k ^ k'. Here 'p"^ is the probability, 
conditioned on b, of finding k photons in the output en- 
vironment modes and is given by 



',(b)\ 



VV = E l^n<ll^ - E 
n:n>k n:n>k 

M r 



E^-n 

n:n>k m=l 



-km)rpk„ 



Again by the NDS state assumption, we have 
(V^f |Vi'^)=/k-<5k,k', 

for 

/k = E |Cn|^(<l)*<l 



n:n>k 



M 



E^-n 

ri:n>k m— 1 



(rori)"'"-'="(ioii)' 



(27) 
(28) 

(29) 
(30) 

(31) 



The orthogonality relations ( 26 ) and ( 29 1 are key to the 



mathematical tractability of the optimal error probability 
problem for NDS input states. 
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A. Optimal (Helstrom) Detection: Error 
Probability and Measurement Operators 

Tracing over the environment modes in (24), we obtain 
Po and pi, which are already in diagonal form by virtue 
of 1261):- 



Ph 



E 



V'k^'' } I have via ( 26 ) and ( 29 1 that 



(32) 



ly ort 



logonal. 



For "Hk = span{V'k°\ 

the 2-dimensional spaces {"Hk} are mutua . 
The joint signal-idler Hilbert space T-l may then be ex- 
pressed as an orthogonal direct sum of the {T^k} with an 
additional component 'H± orthogonal to all the Hk'- 



(33) 



For the purpose of distinguishing po and pi, we may re- 
strict the domain of definition of the difference density 
operator Ap = po — pi appearing in the Helstrom for- 
mula ([5| to 0^ Hk since the support of both po and pi 
is orthogonal to T-L^. Corresponding to (33), Ap may be 



decomposed into a direct sum of operators on Hk 
Ap = 0Ap|j^, Ap|^e/:(Hk), 



where 



(34) 



(35) 



Performing a Gram-Schmidt orthonormalization on each 
•Hk = span{V'^°\?/.^^^} with Vk V II II as the first 

basis vector, we may write the 2x2 matrix of Ap|j^ in 
the Gram-Schmidt basis as [^:- 



Ap 



/ 



(0) 
Pic 



(0) 

Pk 



- I4IVpL"^ 

(0) (1) |7 
Pk Vk - Ml 



Ik 

(0) 
Pic 



l^klVPk 



(0) 



(36) 



On calculating the eigenvalues and trace norm of the 
above matrix, the minimum error probability follows as:- 



^^4-J||e(^Hk) 

k 

k 



(1)\2 
Pk ) 



4|/kP 



1/2 



(37) 
(38) 
(39) 



optimally discriminating po and pi. The optimal mea- 
surement consists of two orthogonal projection operators 
IIo and Hi augmented with an additional projection op- 
erator IIj^ onto "Hj^ to make a complete projective mea- 
surement with Ho + Hi -|- n_L = lis, the identity on "H. 
Note that the state never projects onto 'H±. Hq and Hi 
are given in the usual way [1] by 



and 



no = ^n+{Ap|j 



ni = ^n_{Ap|j, 



(40) 



(41) 



where n_|_/_{Ap|j^} is the projection operator onto the 
one-dimensional eigenspace corresponding to the posi- 
tive/negative eigenvalue of Ap|j^. This space is the span 
of the eigenvector of Ap|j^ of (36 1 with positive/negative 



eigenvalue. In somewhat more physical language, we may 
view the optimum measurement as a two-stage measure- 
ment. In the first stage, we perform a Quantum Non- 
Demolition (QND) measurement, i.e., a measurement 
that projects the state on to one of several orthogonal 
Hilbert spaces without destroying it. In our case, the re- 
ceived state is projected into one of the orthogonal spaces 
■Hk- Depending on the value of k, we further make on "Hk 
the binary Helstrom measurement corresponding to the 
optimum discrimination of IV'k'') and IV^') with condi- 
tional probabihties Pk'V(Pk'^ +Pk ^) and p^^V (Pk''' +Pk ^) 
respectively. From this description, we may expect that 
this measurement is hard to realize in the laboratory 
since the eigenvectors of Ap|j^ are (M' + M)-mode en- 
tangled states in general. 

The case of unequal a priori probabilities, say ttq and 
TTi for Po and pi, may be handled similarly. is now 
given by 



P. 



1, 



FoPo - TTipilli. 



(42) 



The operator ttqPq — nipi also has an orthogonal di- 
rect sum decomposition over the {^k} (which are un- 
changed). However, the matrices for ApL are different 



from (36) and include the prior probabilities. Diagonal- 



izing these matrices yields the optimal error probability 
and the measurement operators. 



Although the above sum may be hard to evaluate ana- 
lytically, numerical computation to any desired accuracy 
is always possible. 

The discussion above also yields the abstract math- 
ematical description of the measurement operators for 



B. Fidelity and Chernoff-type Bounds 



The optimum error probability P^ is usefully bounded by 



the fidelity and Chernoff-type bounds (|6| and 
NDS states, it is straightforward to show, using 



It, 

m 



For 
and 
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( 29 1 in the definition Q , that the fidehty is given by 
•^=(El^k|j (43) 

, k n:n>k m— 1 ^ ^ 

(44) 

The presence of the absolute 
value sign prevents simp ification of the above expression 
without further assumptions, but we note that, for A 
0, (44) is greater and hence more informative, than the 



where we have used (31 



general lower bound (15). 

The Chernoff-type quantities may be likewise com- 
puted to equal 



(45) 



for s £ [0,1]. The 'singularity' in the terms for which 

Pk'' = is only apparent as /k is also zero for those 
terms. Therefore, we need sum only over terms for 
which both p'j^^ ^ 0. We see that Q(0) and Q(l) do 

not equal 1, since |/kp < Pk'Vk^ in general (see [22]). 
This is explained by the fact that the support of po, i.e.. 

Indeed, 



span{i/)^°^} ^ spanji/)^^''}, the support of pi. 



since tp^^ 9^ V'k^^ for k, po a-nd pi do not have the 
same support on "Hk, and consequently also on ®k^k- 
As a result, 0(0) = Tr[pSjpi] = Tr[PoPi] ^ Tr[/pi] = 1, 
for Pq the projection operator onto the support of po and 
I the identity operator on ©k'Hk. Similarly p\ = Pi ^ I 
is the projector onto the support of pi and so 0(1) — 
Tr[po^i] ^ 1. 

For a transmitted state j^I') = ®m=i\'^m) that is a 
product of M signal-idler states iV'm), we have the mul- 
tiplicative properties 



M 



(46) 



IV. QUANTUM READING AND TARGET 
DETECTION 

A. Output State Fidelity and Related Error 
Probability Bounds 

^For the quantum reading and target detection scenarios, 
re have A = 6*1 — 6*0 = 0. In Section II, we obtained the 



lower bound ( 19 ) on the output state fidelity in terms of 



the input state's signal photon probability distribution 
{pn}- For an NDS input state with that {p„}, the fidelity 
( 43 1 evaluates to 



(48) 



which is exactly the general fidelity lower hound (17)- 
(19) of Section II with A = 0. Thus, among all input 



states with a given the NDS states with that 

minimize the fidelity. 

Further, since p = 7'ori + t^ti < 1 when at least one 
Tb < 1, we have by the convexity of the function x ^ p^ 
and Jensen's inequality |l23j that, for a given total signal 
energy TV^, 



' OO 

E 



Priirori -I- toil)" 



> 



(rori + tQti) 



12^=0 "P" 



{rori + toti)^^' 

^min(Afs). 



m—l 



(49) 

(50) 

(51) 
(52) 

The inequality ([50| is an equality (at least at integer 
values of Ns) for precisely the states having Ng total 
signal photons with probability one. The only NDS states 
with this property are the multimode Fock states |A^i) ® 

• • -(^INm) with total number of photons J2m=i ^™ — ^s- 
We have thus shown that a (possibly multimode) Fock 
state with Ns photons minimizes the output fidelity at 
given Ns. 

The above argument, employing convexity as it does, 
does not go through when A 7^ 0. Moreover, even for 
A = 0, achieving minimum fidelity is not equivalent to 
achieving minimum Pe as will be evident in the examples 
to follow in Sections IV.C-E. 



and 



Error Probability Lower Bounds 



M 



m—l 



Since J^min(-^s) of (521 is a lower bound on the output 

(47) 

state fidelity for any multimode pure input state lip) of 
energy Ng, we obtain the universal lower bound 



which simplify computations by converting the sum over 
vector k to a product of scalar sums. We illustrate the 
results of this subsection in Section IV. D by applying 
them to quantmn reading with the EPR state. 



2 (1 - v/l--^min(iV,)j = 



l-(rori+toti)2^= <Pe[|^)] (53) 
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on the output error probabihty Fe[|'/')] using 

For any input state |^') consisting of N copies of an 
(M + M')-mode signal-idler state |?/') 



N 



i*> = (8)i^), 



(54) 



we may also obtain an upper bound on the quantum 



Chernoff exponent ^qcb 



iogQ[^] 

N 



18 . It was shown 



in [TS] that while the finite- behavior of -Pe['3/] is com- 
plicated, Pe [vP] is asymptotically exponential in N in the 
sense that 



(55) 



Assuming has finite average energy A^s (so that \'^) 
has average energy A^ A^s ) , we have from the fidelity lower 
bound J-"min(-^s) < •^['0] coupled with the lower bound of 
^ that 

\ [J'n,UNs)f = |(rori + toii)2(^=)-^ < Pe[f ], (56) 
from which the bound 

Cqcb < log [(rori+ toil) "'^"] (57) 
follows on taking logarithms. 

2. Reading of an Ideal Memory 

For the case of i?i = 1 (called "ideal memory" in [3]), 
the channel £i is simply the identity channel. Therefore, 
the return state pi is pure - to wit, the transmitted state 
In such a situation, the fidelity ./^[l-i/))] equals the 
Chernoff bound Q{\ip)) [T5] , so that we have the bound 
(which is stronger than the upper bound of m)'- 



(58) 



In conjunction with the result of Section IV. A that the 
multimode Fock states attain the minimum fidelity, we 
have the remarkable consequence that among all the 2M- 
mode signal-idler states with signal energy A^^, a signal 
mode Fock state with A^^ total photons has the best (low- 
est) Chernoff bound. See also the discussion in Section 
IV.C.l. 

We mention in this connection that an optimization of 
the Chernoff bound was carried out in [21] fixing the to- 
tal energy in signal and idler but restricting to the cases 
when the input state is a single-mode squeezed thermal 
state and an M = 1 signal-idler two-mode squeezed ther- 
mal state ~ it was shown that the single-mode and two- 
mode squeezed vacuum states minimize the ideal memory 
Chernoff bound in that class. 



3. Target Detection: No- Go Result for Large Loss 

The case of target detection corresponds to Tq = 
and to = 1. Under these conditions, the universal lower 



bound (531 reads 



(59) 



It is easy to check - see (61 ) - that the error probability 



obtained from a pure coherent state of energy A^^ is 



Pe[CS] 



(60) 



The case of large loss i?i ^ 1 is of practical importance 
for standoff target detection. Because 



--Ri 



1 - R 



1 m 



(60) under such conditions, we see that the error proba- 
bility of a general input state of energy N^, which is lower 
bounded by the LHS of ( 59 1 , is not appreciably smaller 
than the coherent state error probability (60 1. We thus 



have a no-go result for appreciable quantum advantage 
in target detection under high loss conditions that ap- 
plies to any multimode input state. Note that this does 
not contradict the 6 dB advantage in the error exponent 
of the EPR state over coherent states claimed in [2] for 
high loss target detection because that analysis was car- 
ried out with the additional assumption of large thermal 
background noise in each signal mode. 

The question arises if one can connect the fidelity with 
the Chernoff bound as we did in Section IV. A. 2 above. 
If the input state has signal-idler entanglement, _Ro = 
and i?i 7^ 1 imply that both po and pi are mixed in gen- 
eral so that the argument connecting the fidelity and the 
Chernoff bound no longer applies. However, for the case 
of a pure transmitted state that is not entangled to any 
idler modes kept at the receiver - this is called a Type 
I target detection scenario in |25| - po is a pure state, 
namely the vacuum state of the signal modes. Thus, it 
is again true that in a Type I target detection scenario, 
the number state transmitter yields the lowest Chernoff 
bound among all pure transmitted states of energy A^^. 
When i?i ^ I prevails, the no-go result given above is 
in force even in this case. In such regimes, the coherent 
state performance is little different from the number state 
performance (see (70) below). This holds even for mod- 
erately large Ri, as we will see in Section IV.E.l (Refer 
Fig. 3). 

In the remaining subsections, we compare in detail the 
performance of coherent states, number states, and the 
EPR states of [3]. 



B. Coherent States 

For a coherent state input of energy Ng, say the single- 
mode state (see [3]) of mean amplitude \/T^, the 
optimal error probability ([s]) for discriminating po and pi 
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evaluates to 



Pe[CS] = 2 " Vl - e-('-i-'-o)'^= 



(61) 



The RHS is exactly the lower bound on the optimal error 
probability of any classical state derived in [5] (Theorem 
1 therein). Therefore, pure coherent states are the opti- 
mal classical states for quantum reading in the absence 
of added thermal noise. A physical realization of the op- 
timal error probability (61) is provided by the so-called 



Dolinar receiver 126' and some other receivers exist which 
approximate that performance. We discuss these issues 
briefly in Sec. IV.E.2. 



U 




Number States 



Let us now consider the obtained by transmitting 
number states. For a transmitted Fock state |N) = 
\Ni) (g) • • • (g) |A^„i) with Ns — X]m=i ^'K' it ^e&Ti that 



probabilities:- 




r (b) 

tor pjj. 



a product of binomial 



M 

n 



h 

r\,fYi 



(62) 



The output states pa and pi commute, so that the tech- 
niques of Section III are not required to evaluate the per- 
formance. The optimal quantum measurement is photon 
counting on the individual modes followed by classical 
processing of the count results. It is easy to verify that 
the optimal decision rule is given by 



R'^ ■ T 



say £i 

> 

< 
say £o 



(63) 



for T the observed total photon count (Recall that we 
are assuming _Ro < i?i). This is equivalent to the rule 



FIG. 2: (Color online) Number state vs. Coherent state Cher- 
noff exponent Gain G of eq. ( 69 1 for values of Ro and Ri 
above 0.4. Appreciable gain is seen for larger values of the 
reflectances. 



the distribution of the photons as well. This feature has 
practical implications that we discuss in Section IV.E.2. 
From (62|-(63), Pe is easily computed numerically, as we 



will do in Section IV.E.l. 

We also have the essentially classical Chernoff Bound 



Pe<l[R'o'Rl-''+TtTl- 



(66) 



where 



In 



-Ri ln(i?o/i?i 
TiMTo/Ti) 



In 



TqRi 



TiRq 



(67) 



so that ^iMum is the number state Chernoff exponent. This 
bound is obtained by computing Q{s) directly from the 
definition ([t]) and finding the optimum exponent analyti- 
cally. The above results apply when Rq ^ and i?i ^ 1. 
For a coherent state transmitter, the Chernoff bound 



T 



say £i 
> 
< 
say 



(64) 



where the threshold equals 



(65) 



Note that the decision rule is independent of the actual 
Fock state chosen as long as Ns is fixed. Moreover, we 
can show from ( 62 ) that the probability of counting T 



photons is also independent of the distribution of input 
photons among the modes as long as the total number 
is Ns- Consequently, is independent of the details of 



(68) 



for ^cs the coherent state Chernoff exponent. A use- 
ful measure for quantifying the improvement obtainable 
from a number state transmitter from coherent state per- 
formance is the ratio of their Chernoff exponents, which 
we call the "gain" G - 



G 



1^ 



(69) 



In Fig. 2, we plot G against the reflectances i?o and 
which shows that the gain is appreciably greater than 
unity only when both i?o and i?i are fairly large. Fur- 
ther performance comparisons at varying values of Ng are 
made in Section IV. E. 
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1. Ideal memory and Target Detection 

For the case Rq = (target detection), we have the exact 
result 



^ 1, 



P. 



(70) 



The optimum decision rule in this limit is to declare the 
target present if and only if the total count T > 0. In 
the opposite limit Ri — 1 (the ideal memory of [3]), we 
likewise obtain 



(71) 



and the optimum decision rule is to declare b = 1 if 
and only if T = Ns- The ideal memory case is excep- 
tional in that, irrespective of the particular (pure) state 
transmitted, pi is a pure state - the transmitted state it- 
self. In such a situation, as mentioned in Section IV. A. 2, 
the fidelity equals the Chcrnoff bound [18 . Thus, the 
Chernoff bound reads 



- 2 



(72) 



1 to ( 71 1 , is in fact 



which, on comparing ( 52 ) with ri 
an equality for Fock states. 



D. EPR State 

We now consider the M-mode EPR state with per-mode 
energy iV = A^^/M of [3], which will afford an illustration 
of the techniques developed in Section III. The EPR state 
is an M-fold tensor product of a two-mode squeezed vac- 
uum state which, for a single mode-pair, has the photon 
number representation 



/EPR 



1 



N- 



N 



N + 1 



i/2 



\n)i\n)s. (73) 



|'0)epr is evidently an NDS state. For the case of no 
added thermal noise being considered in this paper, it was 
shown in [3] using the Bhattacharyya bound (5(1/2) that, 
for Ns > iVth(i?o: -Ri) - a threshold total energy depend- 
ing on the reflectances, there exists an M for which the 
classical error probability is worse than the EPR Bhat- 
tacharyya bound (Theorem 2 of [3]). For the case of ideal 
memory, a similar threshold theorem was shown to hold 
with A^th = 1/2 for all M greater than a threshold M 
(Theorem 3 of [3]). We show here that the techniques of 
Section III may be used to derive these results in an al- 
ternative manner, as well as providing numerical results 
at chosen A^ and M. 

We first obtain the output fidelity and Q{s) for the 
state |i/')epR) i-e., for M = 1. The photon probability 



distribution p^. of ( 27 ) in the output environment mode 



is in the Bose-Einstein form of a thermal state with T^N 
average photons:- 



(b) 

Pk = 



1 



TbA^ + 1 VTbA + l 



Evaluating (30 1 yields 



Ik = 



[(l-rori)A-Hl]^ 



k+l ■ 



We may now use (|43|) and (45) to get:- 

1 



[(l-rori-to^i)A^+l] 



and 



with 



Q{s) 



M=l 



a = > 1, 

TiN + 1 



C = 



[(l-rori)A^ + l]^ 



and 



D = TiN. 
We may verify that the inequalities 

a < fi 

and 

C>D 



(74) 
(75) 

(76) 

(77) 

(78) 
(79) 
(80) 

(81) 
(82) 
(83) 



hold pTj. For a given TV, obtaining the Chernoff bound 
entails finding the exponent s* that minimizes (77). We 



know that Q{s) is a convex (and continuous [5S]) function 
of s in [0,1] [TSl. It is, from ( (TT] ), evidently also twice 
differentiable in s. Accordingly, two cases logically arise 
depending on the sign of Q'(0):- 

1. If Q'(0) > 0, Q[s) is an increasing function of s, so 
that = 0. 

2. If Q'(0) < 0, then may be found by setting 
Q'{ssf) — 0. If this equation has no solution in 
[0,1], = 1. 



Differentiating (77), we find that 

Q'{s) = {Ca' - DI3''Y\\nP°)P'' - (Ina'^)a"]. (84) 
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Consequently, the condition for deciding among the 
above two cases is 



Case 1 



A-^"" ^ 1 



(85) 



Case 2 



In the event of Case 2, setting Q'(s^) = gives as 



In 




In 


a 





(86) 



In the event that the RHS of (86) is greater than 1, we 
have s* = 1. 

Finally, from the multiplicative properties of the fi- 
delity and the Chernoff-type quantities, we have that, for 
the M-mode input state lip) = <8)m=ilV')EPR, the fidelity 
equals 



-2M 



(87) 



^[EPR] = [il-ron-toh)N+l] 
and the Chernoff bound is 

for the s* obtained from the above case analysis. 

When the input state \ip) is a multimode Gaussian 
state, i.e., a state whose Wigner function is a Gaussian 
probability distribution (see e.g. ,[29]), so are po and pi 

since £^ (8) id is a linear, and hence Gaussian, chan- 
nel. The input state |V')epr is a Gaussian state. We 
could therefore use the general Gaussian state technique 
of (30j to derive the fidelity between po a-nd pi . Similarly, 
the technique of symplectic diagonalization [19, may be 
used to derive Q{s), as it was in [3 . 

Let us now connect these results to Theorems 2 and 3 
of |3] - the statements of these theorems were reviewed in 
the first paragraph of this subsection. Theorem 2 makes 
essential use of the large M limit of the Bhattacharyya 
bound (5(1/2) given by ([7|. Some lengthy but straight- 
forward algebra verifies that the M — >■ oo limit (at con- 
stant Ns — MN) of the A/-pair Bhattacharyya bound 



B{N, 



- lim 

2 M-j-oc 



[Q(i/2)Uf. 



(89) 



for Q{1/2)\m=i obtained from (77) is identical to the cor- 
responding result (117) of [31] after adjusting for differing 
notation (see [35]). The same reasoning in ^1] that es- 
tablishes the "Threshold energy" theorem (Theorem 2 of 
[3]) may then be carried out from ( 89 ) to give an identical 
threshold energy as before. 

Consider now the case of an ideal memory with Ti — 
for comparison to Theorem 3 of [3]. That theorem de- 
pends on the expression (129) of [31] for the Chernoff 
bound which we rederive using the method of this sec- 
tion. The same analysis as in [3T can then be used to 



reproduce the result of Theorem 3. From (81), we have 
D = so that ( pTTj ) is clearly minimized at = 1 so that 
the M-mode Chernoff bound is 

'Pe < ^Q[idcal] = i[(l-ro)7V+l]-2A^ (90) 

which agrees with (129) of [ni]-[Sl]- Further, we see that 
the fidelity (87) also equals Q [ideal] as it should because 



Pi is a pure state. In the general case of i?i ^ 1, we can 
study the behavior of the Chernoff bound as a function 
of M and Ns by numerically obtaining for each value 
of these parameters. 



E. Comparison of Coherent State, Number State, 
and EPR State Transmitters 

1. Numerical Comparison of Error Probability 

In this subsection, we compare quantitatively the error 
probability performance of the three types of states con- 
sidered in the previous subsections for target detection, 
and reading of non-ideal and ideal memories. The rep- 
resentative plots below show the error probability on the 
y-axis in logarithmic scale against the total average sig- 
nal energy on the x-axis. We assume the number of 
modes M = 50. We reiterate that the number of modes 
has no effect on either the coherent state or number state 
performance, which depend on Ng alone. For the EPR 
state, varying M changes the performance (as given by 
the Chernoff bound), although the change is not appre- 
ciable once M is around 20 — 30. Thus, the plots given 
here are fairly representative of the best possible EPR 
state performance. 

We summarize how the plots were made. For each Ns, 
the universal lower bound of ( 53 ) was plotted in Figs. 2-5. 



The coherent state error probability is given by the closed 
form expression (61 ). The number state error probability 



was calculated numerically for each value of Ns using the 
count probability distribution ( 62 1 and the decision rule 
(63) for the cases of reading of non-ideal memories. The 



number state Chernoff bound (66)-(67) was also plotted 
for these cases. Closed-form expressions (70) and (71) 
were used for the detection and ideal memory cases. The 
EPR state fidelity expression (87 1 in conjunction with ^ 
gives the EPR fidelity lower bound 



hi- Vl--^[EPR]) < Pe[EPR], 



(91) 



Finally, for each value of Ng, the EPR Chernoff bound 
( 88 ) was calculated for the detection and non-ideal mem- 
ory cases using the procedure for obtaining s* described 
in Sec. IV.D. 

In all the figures, we see that the number state and 
EPR state transmitters eventually outperform the coher- 
ent state transmitter, in tune with the conclusions of [3]. 
In the target detection case (Fig. 3), the number state 
also outperforms the EPR state lower bound, although 
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FIG. 3: (Color online) Error Probability Bounds versus Ns 
for target detection with Ro = and Ri = 0.3. The number 
of modes Af = 50 for the EPR state curves. 
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FIG. 4: (Color online) Error Probability Bounds versus Ns for 
quantum reading with Ro = 0.3 and -Ri = 0.6. The number 
of modes M = 50 for the EPR state curves. 

the performance difference between the three states is not 
appreciable. We see that, even for an i?i that is much 
larger than that expected in a realistic target detection 
scenario, the coherent state performance is not apprecia- 
bly worse than the number state performance. As i?o in- 
creases, the perfomance gain over classical increases also 
as evinced in Figs. 4-7. Fig. 4 shows a case where the 
difference Ri — Rq is small and the reflectances them- 
selves are not very high. In such cases, the number state 
and coherent state performances are not appreciably dif- 
ferent. However, the EPR Chernoff bound drops below 
the number state performance for A^^ greater than about 
40 photons. Fig. 5 also shows a case of small i?i — Ro, 
but the reflectances themselves are appreciable. We see 
that the nonclassical transmitters' gain over the coher- 




Number state Chernoff UB 



10-^ ' ' ' ' ' ' ' 1 

10 20 30 40 50 60 70 80 
Total Signal Energy 

FIG. 5: (Color online) Error Probability Bounds versus Na for 
quantum reading with Ro — 0.5 and Ri = 0.75. The number 
of modes M = 50 for the EPR state curves. 
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FIG. 6: (Color online) Error Probability Bounds versus Ns for 
quantum reading with Ro = 0.2 and _Ri = 0.8. The number 
of modes M = 50 for the EPR state curves. 

ent state increases. Further, the crossover between the 
number state performance and the EPR Chernoff bound 
occurs later (at about 65 photons) and the slopes of these 
two curves are less different than in Fig. 4. Fig. 6 repre- 
sents distinguishing channels with large Ri— Rq. We see 
that the nonclassical gain over the coherent state trans- 
mitter is even greater, with the EPR state doing better 
than the number state - the crossover of the number state 
performance with the EPR Chernoff bound now occurs 
at about 20 photons. Finally, Fig. 7 represents the read- 
ing of an ideal memory with Rq — 0.5. The nonclassical 
gain is now very large, and the number state transmitter 
performs better than the EPR state, as evinced by the 
fact that it lies below the EPR state lower bound. All the 
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FIG. 7: (Color online) Error Probability Bounds versus Ns 
for reading of an ideal memory with Ro — 0.5 and _Ri = 1. 
The number of modes M — 50 for the EPR state curves. 



plots are consistent with the number state vs. coherent 
state ChernofF exponent gain of Fig. 2 and Fig. 7 confirms 
and strengthens the conclusion of Section IV. A. 2 that the 
number state transmitter has the lowest ChernofF bound 
for reading of an ideal memory. 



2. Experimental Considerations 

In this section, we indicate some technological consider- 
ations regarding the availability of sources and detectors 
that bear upon the implementation of quantum reading 
with coherent, number, and EPR states. 

While coherent sources of any energy are readily avail- 
able, the optimal Helstrom detector of [2^ is not easy to 
realize as it employs feedback in addition to linear op- 
tics and photodetection. Nevertheless, the Kennedy re- 
ceiver f33' and the so-called Optimum Displacement Re- 
ceiver (ODR) ^34j achieve the same error exponent and 
are more easily implemented since they do not involve 
feedback. Indeed, the ODR was recently demonstrated 
[35] with an overall detection efficiency of ~ 90%. 

A single-mode number state with A^s > 1 is hard to 
generate with current technology. However, we saw in 
Section IV. C that a state consisting of Ng spatial or tem- 
poral modes each of which is in a single photon state, i.e., 
a state of the form ■ -^ll) with A^ total photons, has 
identical performance. As a large variety of single pho- 
ton sources are available (see, e.g., [33]), the generation of 
multimode Fock state with many modes each in a single 
photon state does not appear problematic with current 
technology. The optimal Helstrom measurement is pho- 
ton counting on the individual modes followed by clas- 
sical processing. While high quantum efficiency photon 
number-resolving detectors are still a developing technol- 
ogy (see, e.g., [37]) and require cooling to superconduct- 
ing temperatures, the source described above consisting 



of multiple single-photon states requires only standard 
single-photon avalanche diode (SPAD) [3H] technology 
with either one or many detectors depending on the mode 
implementation. Detector quantum efficiencies and other 
system losses act as multiplicative factors to Rq and Ri 
and do not essentially change the analysis described in 
Section IV. C, as also is the case for the coherent state 
transmitter result of Section IV. B. SPADs, for instance, 
can attain quantum efficiencies of around 77 ^ 0.75. After 
adjusting Rq and i?i to account for system inefficiencies, 
it is fair to say that the optimum Helstrom receiver for 
quantum reading with number states is realizable with 
current technology. 

Generating multimode EPR states of high total energy 
Ns and large M also does not seem to present a huge 
experimental problem (see, e.g., |39j). Unfortunately, 
the optimum Helstrom detector, while given abstractly 



mod 



by (40)- (41 1, is an entangling measurement over the M 

A realiz- 



es and has no known concrete realization, 
able suboptimal measurement involving homodyne detec- 
tion was proposed in [21 [31] and was shown, significantly, 
to also outperform the coherent state transmitter. How- 
ever, it appears that this comparison was made in j3T] 
only for rather high Rq and Ri . The effects of non- unity 
homodyne detector quantum efficiency, analogously to 
the Fock state and coherent state cases above, have also 
not yet been considered in that measurement. Moreover, 
given the comparisons made in this paper, it is of inter- 
est to see how the suboptimal measurement compares to 
the number state performance for reasonable values of 
Ns and M. 



V. CONCLUSION 

The problem of distinguishing two optical beam split- 
ter channels using multimode signal-idler entangled pure 
states was considered. A general lower bound on the 
output state fidelity and minimum error probability for 
any such input was derived. For Number-Diagonal-Signal 
States, series formulae for the optimum error probabil- 
ity, the output state fidelity, and the Chcrnoff-type up- 
per bounds were derived. For quantum reading and tar- 
get detection, for a given signal photon probability mass 
function, the fidelity bound was shown to be attained 
by NDS states, with multimode Fock states minimizing 
the bound for a given total photon number. For reading 
of an ideal memory with arbitrary states and for Type 
I (signal-only) target detection, the number state was 
shown to yield the best Chernoff bound among all states 
of given energy. For target detection under high loss con- 
ditions, a general no-go result for quantum advantage 
over coherent states was obtained. The above results 
were applied to quantitatively studying the performance 
gains over classical states obtainable by number state and 
EPR state transmitters, which were found to outperform 
the classical transmitters to varying degrees over a wide 
range of reflectances. The experimental outlook on re- 
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alizing the optimal measurement for the number state 
transmitter was argued to be good. It is of interest to 
compare the performance, taking into account reahstic 
experimental parameters, of the non-Helstrom measure- 
ment on the EPR state transmitter suggested in |3] and 
the number state transmitter suggested here. Finally, the 
techniques developed here are likely to prove useful for 
other interesting problems fitting the same framework, 
e.g., the lossy discrimination of optical phase shift chan- 
nels. 
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